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ABSTRACT 


Processing  network  problems  are  minimum  cost  network  flow  problems  In 
which  the  flow  entering  (or  leaving)  a  node  may  be  constrained  to  do  so  In 
given  proportions.  These  problems  Include  many  practical  large-scale  linear 
programming  applications.  In  this  paper  a  special  class  of  processing 
network  problems  called  allocation/processing  (AP)  network  problems  Is 
described.  AP  network  problems  model  a  realistic  situation  In  which  a 
single  raw  material  is  allocated  to  factories  and  the  resulting  finished 
products  are  distributed  to  customers. 

A  special  partitioning  method  (PPR  algorithm)  is  introduced  for  the 
solution  of  AP  network  problems.  This  algorithm  maintains  a  working  basis 
as  well  as  a  so-called  representative  spanning  tree.  A  method  for  updating 
the  representative  spanning  tree  is  discussed  which  avoids  tracing  a  number 
of  cycles  at  each  iteration. 

A  FORTRAN  Implementation  (PPRNET)  of  the  PPR  algorithm  Is  described. 
Comparative  computational  results  are  presented  for  randomly  generated  AP 
network  problems  which  were  solved  by  PPRNET  and  MINOS.  These  results 
indicate  that  significant  computational  gains  are  possible  with  the  use  of 
special  purpose  processing  network  codes. 


Key  Words:  Network  Flow  Algorithms,  Linear  Programming,  Large  Scale 


Systems,  Processing  Networks. 
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1 .  INTRODUCTION 

Processing  network  problems  are  a  special  type  of  embedded  network 
problem.  For  these  problems,  the  network  has  two  special  types  of  nodes 
called  refining  nodes  and  blending  nodes.  At  a  refining  node,  the  entering 
flow  splits  Into  several  leaving  components  In  fixed  proportions  while  at  a 
blending  node,  several  components  enter  In  fixed  proportions  to  form  the 
leaving  flow.  The  term  processing  node  is  used  to  describe  either  a 
refining  node  or  a  blending  node.  A  number  of  applications  which  exhibit 
processing  network  structure  have  been  listed  by  Roene  [ 18] ,  including 
energy  modelling  [20]  and  assembly  models  [25].  McBride  [21]  notes  that 
short-term  financial  planning  models  [19]  can  be  formulated  as  processing 
network  problems.  Charnes  et  al.  [5]  have  recently  described  a  processing 
network  model  which  Is  used  in  manpower  planning. 

The  computational  study  [11]  showed  that  special  purpose  network  codes 
are  150-200  times  faster  than  the  APRX  III  linear  programming  (LP)  code. 
This  raises  the  question  of  whether  It  Is  possible  to  develop  special 
purpose  processing  network  codes  which  are  significantly  faster  than 
state-of-the-art  LP  codes.  We  believe  the  answer  to  this  question  is  in  the 


affirmative,  and  we  present  computational  results  In  this  paper  which 
support  this  position. 

Although  processing  network  problems  can  be  solved  by  embedded  network 
techniques  as  described  by  Glover  and  Klingman  [12]  and  Gupta  and  McBride 
[13],  there  have  been  specialized  algorithms  proposed  by  Koene  [18]  and 
McBride  [21]  which  take  advantage  of  processing  network  structure. 
McBride's  algorithm  has  been  implemented,  and  some  preliminary  computational 
results  are  presented  in  [21].  Charnes  et  al.  [5]  describe  a  heuristic 
procedure  for  the  solution  of  processing  networks.  Beck,  Lasdon  and 
Engqulst  [4]  solve  network  problems  with  equal  flow  constraints  by  means  of 
a  penalty  method.  In  this  paper  we  describe  yet  another  algorithm  for 


solving  certain  processing  network  problems,  and  we  discuss  its 
implementation  and  testing. 


2 .  TEST  PROBLEMS 

All  processing  nodes  for  the  class  of  problems  for  which  computational 
tests  were  conducted  are  refining  nodes.  For  processing  node  j,  as  shown  in 
Figure  1,  the  flow  on  an  arc  (j,k)  must  be  a  fraction  of  the  flow 

entering  node  j.  Since  conservation  of  flow  must  hold  at  node  j,  it  is 
required  that 

P 

I  ° jk(v>  -  1  (2.1) 

v«l 

We  note  that  processing  nodes  are  represented  graphically  as  squares  while 
ordinary  nodes  are  shown  as  circles.  In  Figure  1,  we  refer  to  the  arcs 
leaving  node  j  as  processing  arcs  and  to  the  arc  entering  node  j  as  an 
allocation  arc. 


Figure  1.  A  Processing  Node. 
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Processing  network  problems  can  be  formulated  as  networks  with  either 
side  constraints  or  side  variables,  and  we  proceed  to  describe  these  two 
tormulatlons  as  they  apply  to  our  class  of  test  problems.  The  first  of 
these,  called  the  sparse  formulation,  is  used  by  our  solution  algorithm 
which  is  described  in  Section  4.  In  addition  to  the  usual  conservation  of 
|  flow  constraint  at  processing  node  j  of  Figure  1,  the  sparse  formulation 

requires 

j  ^jk(v)  xjk(v)  ”  ®Jk(v+l)  xjk(v+l)  “  v»l ,2, . . . ,p-l  (2.2) 

where  -  ajk(v)  and  xJk(v)  is  the  flow  on  arc 

The  second  formulation  is  known  as  the  compact  formulation.  In  this 
formulation,  the  processing  arcs  of  Figure  1  are  replaced  by  a  single 
processing  column  in  the  constraint  matrix  of  the  form 

1  in  row  j 

~a jk(v)  in  row  v-1»2 . .  (2.3) 

|  0  elsewhere . 

In  the  computational  tests  described  in  Section  3,  the  compact  formulation 
Is  used  to  generate  an  input  file  for  the  NINOS  optimization  code  of  Murtagh 
and  Saunders  [22]. 

A  further  requirement  for  our  test  problems  is  that  all  allocation  arcs 
emanate  from  a  common  node.  These  test  problems  are  called 
al location/ processing  (AP)  network  problems,  and  they  may  be  described  as 
follows:  A  certain  amount  (s)  of  raw  material  is  supplied  at  the  allocation 


S 


node.  This  raw  material  Is  to  be  allocated  to  factories  (processing  nodes) 
which  all  produce  the  same  finished  products  from  the  raw  material. 
However,  the  proportion  of  a  unit  of  raw  material  devoted  to  a  particular 
product  varies  from  factory  to  factory.  After  a  finished  product  leaves  the 
factory.  It  enters  a  transportation  subnetwork  devoted  to  distribution  of 
this  type  of  product.  All  such  transportation  subnetworks  are  identical  and 
their  destination  nodes  have  customer  demands  (d^)  which  must  be  satisfied. 

The  topology  of  an  AP  network  is  determined  by  the  following 
parameters:  the  number  of  processing  nodes  r,  the  number  of  finished 
products  p,  and  the  number  of  destinations  q  in  a  transportation  subnetwork. 
An  example  of  an  AP  network  with  r  »  2,  p  *  2,  and  q  *  2  is  shown  in 
Figure  2.  In  this  example,  node  1  is  the  allocation  node,  nodes  2  and  3  are 
the  processing  nodes,  and  nodes  4  through  11  are  contained  in  the 
transportation  subnetworks. 


Figure  2.  An  AP  Network  Example. 
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The  allocation  and  processing  arcs  in  AP  problems  are  capacitated  while 
all  other  arcs  are  uncapacitated.  The  capacities  on  the  allocation  arcs 
model  the  limited  throughput  of  the  factories.  No  capacities  were  Imposed 
on  the  arcs  of  the  transportation  subnetworks  since  the  density  of  these 
subnetworks  indicates  abundant  means  of  transportation. 

It  follows  from  the  theory  of  transportation  problems  that  for  any 
feasible  solution  of  an  AP  network  problem,  the  total  flow  into  a 
transportation  subnetwork  must  be  equal  to  the  total  demand  for  that 
subnetwork.  The  amount  of  flow  into  a  given  transportation  subnetwork  is 
also  equal  to  a  sum  of  terms  of  the  form  x  (flow  on  corresponding  allo¬ 
cation  arc).  This  gives  rise  to  p  equations  in  r  unknowns.  In  problems 
where  the  aj^  are  randomly  generated  and  r  <  p,  this  system  of  equations 
will  have  a  unique  solution  almost  surely  (if  the  dt  are  chosen  to  guarantee 
that  the  system  is  consistent).  If  such  a  unique  solution  exists,  the  AP 
network  problem  can  be  decomposed  into  p  transportation  problems.  For  this 
reason,  we  restricted  computational  tests  to  AP  networks  with  p  <  r. 

3.  BASIS  STRUCTURE 

We  now  consider  the  LP  problem  which  corresponds  to  an  AP  network 


problem  with  parameters  p,  q,  and  r  using  the  sparse  formulation. 


/V 
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,v  v.r. 


i 


minimize 


subject  to  Ax  »  b 


0  <  x  <  u 


(3.1) 


In  (3.1),  A  is  an  mxn  matrix;  c,  x,  and  u  are  n-vectors;  and  b  is  an 
m-vector.  If  an  arc  is  to  be  uncapacitated  this  is  accomplished  by  choosing 
a  sufficiently  large  component  of  u  for  that  arc.  Without  loss  of 
generality,  it  is  assumed  for  each  processing  node  as  shown  in  Pigure  1, 


ljk(v)  uij  M  ujk(v)  »  v  *  1»2*.“»P* 


(3.2) 


The  matrix  A  in  (3.1)  can  be  written 


(3.3) 


If  we  let  raj  *  m  -  r(p  -  1)  and  nj  »  n  -  rp,  then  the  mj  x  nj  matrix  Aj  is 
the  node-arc  Incidence  matrix  for  the  pure  network  (non-processing)  arcs  of 
the  problem.  The  raj  x  rp  matrix  A£  is  the  node-arc  incidence  matrix  for  the 
processing  arcs,  and  the  r(p  -  1)  x  rp  matrix  A3  contains  the  additional 
proportional  flow  constraints  of  (2.2).  We  assume  without  loss  of 
generality  that  an  extra  column  containing  a  single  1  with  all  other 
components  equal  0  (root  arc  column)  has  been  adjoined  to  Aj  and  that  the 
rank  of  Aj  is  mj .  It  is  easily  deduced  that  the  rank  of  A3  is  r(p  -  1)  and 
hence  that  the  rank  of  A  is  m. 
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Froo  now  on,  B  denotes  a  basis  of  (3.1).  Thus,  B  is  a  nonsingular 
m  x  m  submatrix  of  A. 

Each  processing  node  j  as  shown  in  Figure  1  has  p  processing  arcs 
emanating  from  it.  It  follows  from  (2.2)  that  B  contains  either  p  or  p-1  of 
the  columns  corresponding  to  these  processing  arcs.  If,  for  processing  node 
j,  B  contains  p  processing  arcs,  then  processing  node  j  is  said  to  be 
active.  Otherwise,  processing  node  j  is  inactive.  If  a  spanning  tree  of 
the  AP  network  has  the  property  that  it  contains  a  single  processing  arc  for 
each  active  processing  node,  it  is  called  a  representative  spanning  tree 
(RS-tree).  If,  further,  all  allocation  arcs  are  included  in  an  RS-tree,  we 
call  such  a  tree  a  representative  spanning  allocation  tree  (RSA-tree).  In 
either  case,  the  processing  arc  corresponding  to  processing  node  j  is  said 
to  be  the  representative  arc  of  node  j. 


Theorem  1  [18].  For  every  basic  feasible  solution  of  (3.1),  there  exists  a 
corresponding  basis  B  such  that 


B  - 


T 

D 


C 

F 


(3.4) 


where  the  m^  x  m^  matrix  T  is  the  node-arc  incidence  matrix  of  an  RS-tree 
and  the  mj  x  r(p-l)  matrix  C  contains  columns  of  A2. 


Corollary.  In  Theorem  1,  we  may  conclude  that  T  is  the  node-arc  incidence 
matrix  of  an  RSA-tree. 

Proof.  Starting  with  B  from  Theorem  1,  we  can  make  degenerate  pivots 


using  (3.2)  until  a  new  B  is  obtained  for  which  T  in  (3.4)  is  as  desired 


The  following  result  is  useful  whenever  a  partitioned  basis  matrix  is 
maintained.  A  proof  may  be  found  in  [17]. 

Theorem  2.  If  the  matrices  of  (3.4)  are  given  and  the  matrix  Q  is  defined 
by 


Q  -  F  -  DT_1C  , 


(3.5) 


then  Q  is  nonsingular  and 


T-1  +  x-lCQ“'1DT~1  -T-1CQ_l 


-Q-1DT-1 


,-l 


(3.6) 


The  tree  structure  associated  with  T  allows  computations  involving  T-* 
to  be  carried  out  using  highly  efficient  tree  labeling  procedures  [2]  so 
that  T~*  need  not  be  explicitly  maintained.  The  matrix  Q~*  must  also  be 
maintained  in  some  form  in  order  to  carry  out  the  computations  required  by 
the  simplex  method.  However,  Q  is  an  r(p  -  1)  dimensional  matrix,  and  for 
purposes  of  computational  testing,  we  maintained  a  smaller  matrix  R  which 
has  dimension  r.  The  matrix  R  is  constructed  from  Q  by  using  certain  row 
operations  which  utilize  the  special  structure  that  Q  inherits  from  (2.2). 
We  omit  the  details  of  this  construction  but  we  note  that  in  our 
computational  testing  the  matrix  R  was  kept  in  dense  form  along  with  its 
inverse.  The  multiplication  of  a  row  or  column  vector  by  Q""^  is  replaced  by 
a  corresponding  multiplication  involving  R“^.  In  other  words,  for  the 


algorithm  described  in  the  next  section,  R  plays  the  role  of  a  working 


4.  A  SIMPLEX  VARIANT  AND  ITS  IMPLEMENTATION 


The  steps  of  the  revised  simplex  algorithm  are  well  known.  These  steps 
can  be  specialized  to  a  situation  in  which  a  partitioned  basis  is  maintained 
for  a  network  problem  with  side  constraints.  Such  a  specialization  is 
described,  for  example,  in  Chapter  7  of  [17].  In  this  section,  we  point  out 
some  of  the  differences  between  our  algorithmic  approach  and  that  of  [17]. 
Our  algorithm,  which  is  restricted  to  AP  network  problems,  may  be  described 
as  a  partitioning  algorithm  for  certain  processing  network  problems  (PPR 
algorithm).  We  also  discuss  a  FORTRAN  implementation  of  this  algorithm 
called  PPRNET. 

It  will  be  useful  to  make  the  convention  for  AP  networks  that  the 
allocation  node  is  node  1  while  the  processing  nodes  are  nodes  2  through 
r  +  1 . 

In  the  initialization  step  of  the  PPR  algorithm,  we  set  up  an  initial 
RSA-tree  and  the  corresponding  matrix  R.  This  RSA-tree  is  constructed  using 
the  minimum  average  cost  procedure  which  is  described  after  we  introduce 
some  notation.  The  forward  set  F(u)  of  a  node  u  in  a  network  is  defined  to 
be  the  set  of  all  nodes  v  such  that  an  arc  (u,v)  of  the  network  exists.  The 

A 

average  cost  cj  of  processing  node  j  in  an  AP  network  is  defined  by 


:lj  +  ~  l  °  Ik  l  cku 
kcF(j)  ueF(k) 


(4.1) 


It  should  be  noted  that  (4.1)  represents  the  expected  cost  incurred  by  a 
unit  flow  through  node  j  when  the  fraction  of  this  flow  reaching  a 
transportation  subnetwork  is  equally  likely  to  take  any  available  arc.  The 
procedure  for  obtaining  the  minimum  average  cost  RSA-tree  follows. 


*_  * _ N. 
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(i)  Select  processing  node  j'  for  which 

cw  ■  min  Ci.  (4.2) 

J  2<J<r+l  J 

(ii)  Introduce  p  artificial  arcs  into  the  RSA-tree  each  emanating  from  the 
allocation  node  and  terminating  at  a  distinct  node  of  F(j').  The  flow 
on  each  such  arc  equals  the  total  demand  for  the  transportation 
subnetwork  with  which  it  is  incident.  These  arcs  are  uncapacitated 
with  a  Big-M  cost. 

(iii)  For  each  k  e  F(j'),  include  q  arcs  (k,u),  u  e  F(k)  in  the  RSA-tree. 
Arc  (k,u)  has  flow  du. 

(iv)  For  each  origin  node  k  of  a  transportation  subproblem  which  is  not  in 
F( j ' ) ,  choose  an  arc  (k,u')  such  that 

cku'  "  oin  x  cku  *  (4*3) 

ueF(k) 

These  arcs  are  included  in  the  RSA-tree,  with  zero  flow. 

(v)  Include  all  allocation  arcs  in  the  RSA-tree. 

To  accompany  this  RSA-tree,  it  is  necessary  to  specify  for  each  processing 
node  a  processing  arc  which  will  be  nonbasic  (all  processing  nodes  are 
initially  inactive).  These  processing  arcs  can  be  chosen  arbitrarily.  Once 
they  are  chosen,  the  matrix  R,  which  is  initially  diagonal,  is  determined. 

Pricing  by  the  PPR  algorithm  is  done  by  first  pricing  the  processing 
arcs  and  then  the  pure  network  arcs.  The  nonbasic  processing  arcs  are 
included  in  a  candidate  list  which  is  maintained  with  a  threshold  strategy 
[9].  According  to  this  strategy,  on  the  first  pass  of  the  candidate  list, 
the  pivot  eligible  arc  (j'.k')  with  the  largest  absolute  reduced  cost  is 


taken  as  the  entering  arc*  On  subsequent  passes,  an  entering  arc  (j,k)  oust 
satisfy 

lcjki  >  yfcj'fc'l  <4*4> 

where  Cjk  and  cj'k'  denote  reduced  costs  and  the  parameter  y  Is  a  fraction. 
In  PPRNET,  0.5  Is  the  value  of  y .  When  the  candidate  list  yields  no  more 
pivot  eligible  arcs,  the  pure  network  arcs  are  priced  using  the  outward  node 
most  negative  rule  [10],  [24].  After  the  pure  network  arcs  have  been 
priced,  the  processing  arcs  are  again  priced,  and  so  on  until  no  pivot 
eligible  arcs  remain. 

We  note  that  when  the  ratio  test  is  carried  out,  the  allocation  arcs 
are  not  considered  since  they  never  leave  the  basis.  This  feature  of  the 
PPR  algorithm  is  made  possible  by  (3.2). 

The  basis  updating  procedure  for  the  PPR  algorithm  is  quite  different 
from  that  of  the  algorithm  described  In  [17],  The  RSA-tree  structure  must 
be  maintained,  and  there  are  two  cases  to  be  considered.  In  the  first  case, 
the  leaving  arc  lies  on  the  cycle  formed  in  the  RSA-tree  by  the  entering 
arc,  and  no  allocation  arcs  lie  on  this  cycle.  In  this  case,  the  RSA-tree 
is  updated  as  in  an  ordinary  network  pivot  and  the  matrix  R  needs  no 
updating.  In  the  remaining  case,  it  may  be  required  that  several  subtrees 
be  rehung.  We  consider  the  current  set  of  representative  arcs  to  be  a 
matching  between  active  processing  nodes  and  subtrees.  The  set  of  subtrees 
is  changed  when  the  leaving  arc  is  deleted  and  the  entering  arc  is  adjoined. 
If  we  put  a  cost  of  0  on  the  current  representative  arcs  and  a  cost  of  1  on 
all  other  processing  arcs,  then  the  problem  of  updating  the  RSA-tree  with  a 
minimum  number  of  changes  in  representative  arcs  becomes  one  of  finding  a 
minimum  cost  augmenting  path.  Augmenting  path  methods  are  covered  in  [8]. 


In  order  to  obtain  an  augmenting  path,  PPRNET  uses  Dial's  Implementation  [6] 
of  the  shortest  path  algorithm  of  Dljkstra  [7].  Once  the  RSA-tree  has  been 
updated,  the  matrix  R  is  regenerated.  PPRNET  then  computes  R”*  using  the 
IMSL  [16]  subroutine  LINV1F.  The  matrix  R~*  is  handled  in  this  way  since 
PPRNET  is  intended  for  experimental  use  only. 

The  data  structure  of  PPRNET  is  patterned  after  those  of  pure  network 
codes.  The  RSA-tree  has  the  predecessor,  depth,  thread  and  reverse  thread 
[2]  associated  with  it.  Additional  arrays  are  used  to  carry  out  the 
operations  required  in  using  the  partitioned  basis  inverse  (3.6).  PPRNET 
uses  five  arrays  of  length  n;  however,  two  of  these  were  included  for 
convenience  and  could  be  eliminated  without  affecting  the  performance  of  the 
code.  Also  included  are  eleven  arrays  of  length  m^,  eight  arrays  of  length 
pr,  eleven  arrays  of  length  r,  and  the  two  r  x  r  matrices  R  and  R-*.  The 
use  of  network  data  structures  in  processing  network  codes  is  Important 
since  it  allows  extremely  large  problems  to  be  solved  in  main  memory. 

5.  COMPUTATIONAL  TESTING 

A  FORTRAN  code  called  PRGEN  was  developed  which  generates  random  AP 
network  problems.  PRGEN  creates  two  files — corresponding  to  the  two  problem 
formulations — for  a  given  problem.  In  order  to  insure  that  the  problems 
generated  are  feasible,  PRGEN  creates  a  feasible  flow  as  it  builds  the 
network.  The  capacity  on  an  allocation  arc  is  set  to  be  y  times  the 
feasible  flow  generated  for  that  arc.  The  value  y,  which  is  at  least  one, 
represents  the  excess  capacity  of  the  factory  system,  and  it  is  an  input 
parameter  for  PRGEN.  For  a  given  processing  node  j,  all  but  one  are 
chosen  randomly  from  a  certain  problem  dependent  range  with  the  remaining 
value  chosen  so  that  condition  (2.1)  holds. 
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M1H0S  [22],  Is  designed  to  solve  linearly  constrained  problems  with 
either  linear  or  nonlinear  objectives.  However,  we  only  made  use  of  the 
linear  programing  capabilities  of  MINOS.  MINOS  uses  the  revised  simplex 
method  with  Bartels-Golub  updating  [3]  of  a  factored  basis  inverse.  The 
basis  reinversion  routine  Is  derived  from  the  method  of  Hellerman  and  Rarlck 
[14],  [15].  All  parameters  for  MINOS  were  set  to  default  values  [23]  except 
for  the  partial  price  parameter.  This  parameter,  which  we  denote  as  k  in 
Table  2,  is  used  to  partition  the  nonbaslc  variables  into  k  subsets  of  equal 
size  for  purposes  of  pricing.  The  starting  basis  for  MINOS  is  created  by  a 
routine  which  selects  from  all  columns  of  the  constraint  matrix. 

For  PPRNET,  in  all  occurrences  where  testing  for  sero  was  required,  a 
tolerance  of  10“*  was  used.  The  Big-M  value  used  in  PPRNET  was  99999. 
Although  no  numerical  difficulties  were  encountered  with  the  Blg-M  Method, 
It  seems  wise  to  replace  it  In  future  testing  with  a  Phase  I,  Phase  II 
approach  or  with  an  exact  non-Archimedean  approach  such  as  in  the  NONARC 
code  [ 1 ] . 

The  parameter  values  for  the  seven  problem  groups  tested  are  shown  in 
Table  1.  The  number  of  problems  in  each  group  is  also  shown.  Groups  four 
and  five  contain  Identical  problems  except  that  the  capacities  on  the 
allocation  arcs  are  larger  In  group  five.  The  costs  for  the  problems 
generated  are  randomly  chosen  and  for  the  allocation  arcs  they  are  between 
500  and  1000  while  for  the  arcs  In  the  transportation  subnetworks  they  lie 
between  1  and  100.  The  costs  on  processing  arcs  are  zero.  The  total  supply 
for  these  problems  varied  from  500  to  900. 

The  test  problems  were  solved  using  the  MNF  compiler  on  the  CDC  Dual 
Cyber  170/750  at  the  University  of  Texas.  The  results  as  shown  In  Table  2 
are  averages  for  each  problem  group.  The  execution  times  reported  are  In 
central  processor  seconds  and  are  exclusive  of  input  and  output.  Invert 
time  Is  the  total  time  required  by  the  LINV1F  subroutine  in  computing  R“*. 
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The  non-netwdrk  Iterations  are  those  iterations  of 

the  PPR  algorithm  which  1 

'•  •• 

p.V- 

Involve  the  regeneration  and 

relnversion  of  R.  Also  reported 

in  Table  2  is 

i 

the  number 

of  processing  nodes 

active  at  optimality 

It  should  be  noted  that  subsequent  to  the 

generation 

of  the  test 

results  In 

Table  2 , 

we  found 

that  MINOS  times  could  be  decreased  by  about  I 

v"* 

1 

1 

25%  in  problem  groups  one 

through 

five  by  setting  the 

partial  price  1 

parameter 

to  30.  Presumably, 

a  similar  decrease  is  possible  In  group  six  as  1 

1 

TABLE  1 

AP 

NETWORK  PROBLEM  DATA 

Problem 

No. 

Compact  Formulation  f 

Group 

Problems  r 

P 

q  u 

No.  Rows 

No.  Columns 

s 

• 

1 

5 

10 

5 

40  1.1 

261 

2020 

. 

2 

5 

16 

8 

16  1.1 

273 

2080 

>  / 

3 

5 

20 

10 

10  1.1 

321 

2040 

4 

5 

20 

5 

20  1.1 

221 

2040 

5 

5 

20 

5 

20  1.5 

221 

2040 

i 

6 

2 

20 

15 

15  1.1 

546 

4540 

7 

1 

20 

8 

50  1.1 

581 

8040 

TABLE  2 

1 

PERFORMANCE  STATISTICS 

FOR  AP  NETWORK  PROBLEMS 

Problem 

Method 

Total 

Invert 

Total 

Non-network  Active  1 

Group 

Time 

Time 

Iterations 

Iterations 

Proc. Nodes 

* 

l 

MINOS,  k-2 

49.2 

1040 

1 

PPRNET 

14.5 

3.6 

1548 

910 

5 

2 

MINOS,  k-2 

70.6 

1325 

2 

PPRNET 

36.7 

15.4 

2832 

1385 

8.2 

3 

MINOS,  k-2 

77.6 

1328 

3 

PPRNET 

38.5 

18.4 

2381 

997 

10 

« 

4 

MINOS,  k-3 

52.2 

1155 

4 

PPRNET 

31.3 

18.1 

2022 

985 

5.2 

5 

MINOS,  k-3 

53.6 

1141 

5 

PPRNET 

40.2 

23.5 

2506 

1269 

11 

6 

MINOS,  k-7 

256.0 

2665 

6 

PPRNET 

138.3 

52.0 

6653 

2782 

15 

• 

7 

MINOS 

DID 

NOT 

RUN 

r-r 

7 

PPRNET 

276.7 

113.6 

12208 

6128 

8 

*  See  text 

for  reduced  MINOS 

times 

■ 
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well.  When  this  decrease  is  accounted  for,  the  ratio  of  total  MINOS  times 
to  total  PPRNET  times  ranges  from  about  2.5  for  problem  group  one  down  to 
1.0  for  problem  group  five.  Although  several  parameters  are  varied  over 
this  range  of  problems,  it  seems  clear  that  the  number  of  processing  nodes  r 
plays  a  major  role  in  determining  PPRNET  solution  times.  In  examining  these 
solution  times,  one  should  keep  in  mind  that  for  each  non-network  iteration 
it  is  necessary  for  PPRNET  to  regenerate  and  reinvert  the  r  x  r  matrix  R. 
For  more  advanced  processing  network  codes,  the  time  spent  in  maintaining 
the  working  basis  Inverse  should  be  significantly  decreased.  The  number  of 
processing  nodes  active  at  optimality  varies  from  40Z  to  75Z  of  the  existing 
processing  nodes.  This  is  of  interest  since  for  the  algorithms  described  in 
[12], [18], [21]  the  size  of  the  working  basis  depends  on  the  number  of  active 
processing  nodes. 

6.  CONCLUDING  REMARKS 

We  have  presented  encouraging  computational  results  for  a  specialized 
optimization  code  for  solving  certain  processing  network  problems.  These 
results  indicate  that  it  is  possible  to  achieve  significant  computational 
gains  for  problems  in  which  the  non-network  portion  is  relatively  small. 
The  need  for  developing  such  codes  is  manifested  by  the  widespread 
applications  which  are  possible. 
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Processing  network  problems  are  minimum  cost  network  flow  problems  In 
which  the  flowing  entering  (or  leaving)  a  node  may  be  constrained  to  do  so 
in  given  proportions.  These  problems  include  many  practical  large-scale 
linear  programming  applications.  In  this  paper  a  special  class  of  process¬ 
ing  network  problems  called  allocation/processing  (AP)  network  problems  is 
described.  AP  network  problems  model  a  realistic  situation  in  which  a 
single  raw  material  is  allocated  to  factories  and  the  resulting  finished 
__products_are  distributed  to  customers. 
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20.  ABSTRACT  (continued) 

A  special  partitioning  method  (PPR  algorithm)  is  introduced  for  the 
solution  of  AP  network  problems.  This  algorithm  maintains  a  working  basis 
as  well  as  a  so-called  representative  spanning  tree.  A  method  for  updating 
the  representative  spanning  tree  is  discussed  which  avoids  tracing  a  number 
of  cycles  at  each  iteration. 

A  FORTRAN  implementation  (PPRNET)  of  the  PPR  algorithm  is  described. 
Comparative  computational  results  are  presented  for  randomly  generated  AP 
network  problems  which  were  solved  by  PPRNET  and  MINOS.  These  results 
indicate  that  significant  computational  gains  are  possible  with  the  use  of 
special  purpose  processing  network  codes. 
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